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ABSTRACT 

We show how to measure cosmological parameters using observations of 
inspirahng binary neutron star or black hole systems in one or more gravitational 
wave detectors. To illustrate, we focus on the case of fixed mass binary systems 
observed in a single Laser Interferometer Gravitational-wave Observatory 
(LIGO)-like detector. Using realistic detector noise estimates, we characterize 
the rate of detections as a function of a threshold signal-to- noise ratio po, the 
Hubble constant Hq, and the binary "chirp" mass. For po — Hq — 100 
km/s/Mpc, and IAMq neutron star binaries, the sample has a median redshift 
of 0.22. Under the same assumptions but independent of i^o, a conservative rate 
density of coalescing binaries (8 x 10~^ yr~^ Mpc~^) implies LIGO will observe 
~ 50yr~^ binary inspiral events. 

The precision with which Hq and the deceleration parameter qq may be 
determined depends on the number of observed inspirals. For fixed mass binary 
systems, ~ 100 observations with po = 10 in the LIGO detector will give Hq to 
10% in an Einstein-DeSitter cosmology, and 3000 will give go to 20%. For the 
conservative rate density of coalescing binaries, 100 detections with po = 10 will 
require about 4 yrs. 

Subject headings: gravitation, binaries: close, cosmology: theory, distance scale, 
methods: statistical 
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1. Introduction 

When completed in the late 1990's, the Laser Interferometer Gravitational-wave 
Observatory (LIGO) (Abramovici et al. 1992) will consist of two interferometers. 
Concurrently, the VIRGO consortium (Bradaschia et al. 1990) will complete a single 
interferometer of comparable precision. In this letter we introduce a new and general class 
of cosmology tests based on the anticipated observation of the gravitational radiation from 
inspiraling binary neutron star and/or black hole systems. 

Schutz (1986) and Krolak and Schutz (1987) noted that observation of binary inspiral 
in three independent interferometers will reveal the source's luminosity distance d^. With 
an independent measure of the source rcdshift, such observations can determine cosmological 
parameters (e.g., Hubble's constant Hq). Our tests differ from theirs: the simplest can 
determine cosmological parameters with observations made in a single interferometer and 
without any other independent knowledge about individual binary systems. Extensions 
systematically exploit additional information available from multiple interferometers. 

To assess the capabilities of these tests we consider observations in a single LIGO-hke 
Interferometer (LLI); thus, our work applies directly to the individual interferometers in the 
LIGO and VIRGO projects. In addition, the two components making up the LIGO detector 
(LD) will be arranged in nearly the same plane and with nearly identical arm orientations 
(Whitcomb 1992) - a configuration which operates like a single LLI of increased sensitivity 
(Finn and Chernoff 1993, FC). We apply our recent results (FC) for the precision with 
which a binary's parameters can be determined in a single LLI. 

Below we describe the physical basis and methodology of our tests. As a concrete 
application we specialize to neutron star binary systems of constant mass. We find the 
effective volumes sampled and the cosmological source rates. We derive the precision with 
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which Hq and go Kiay be determined in a Friedmann- Robert son- Walker (FRW) universe. 
We demonstrate that source evolution may be handled in a straightforward manner. 



2. Gravitational radiation from inspiraling binaries 

Assume that quadrupole-formula gravitational radiation adiabatically damps a 
binary's Newtonian orbit.0 Let the distance between the source and the LLI be d. The 
time-dependent response to the gravity-wave is (G = c = 1) 

m{t) = ^ {irMff e (ns, n^) cos (^f 271 fdt + ^) , (1) 

where the "chirp" mass M = /i^/^M^/^ and /i and M are the reduced and total mass of the 
binary system. The frequency / = f{Ai, T — t) where T is the instant the binary separation 

given the above approximations and \1/ is a constant phase. B is a known function 
of four angles specifying the orientation of the binary system with respect to the LLI 
(0 < < 4, PC): fiB lies along the binary's orbital angular momentum and "fis points from 
the detector to the binary. The wave frequency / is twice the binary's orbital frequency, 
and a LLI will be sensitive to the radiation during the last few minutes of inspiral as / 
sweeps from ~ 10 Hz to ~ 1 KHz. For binaries at cosmological distances the form of the 
response is the same, only now in place of d we have the luminosity distance di and in place 
of M. we have }A' = (1 + z)A4, where z is the redshift. For more detail, see FC. 

^Cutler et al. 1992 have shown that higher order relativistic contributions modify the 
orbital evolution. The tests described herein are entirely compatible with the use of a more 
accurate evolution, and we do not expect our conclusions to be altered significantly by the 
adoption of a refined waveform. 
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Schutz (1986) showed that three interferometers of different orientations are capable 
of determining rf^, A4', ^, T, fis, and n^. In contrast, a single LLI can measure only 
A = Q/di, M.', T, and "if (FC); thus, (II cannot be determined directly since hs and fiB 
are unknown. However, we know that fis and fiB are each uniformly distributed on the 
sphere; consequently, the distribution of ©(ng, 7x3) is also known. Thus, the observed 
distribution of A is related to the number density of binary systems on spheres of constant 
(II- Similarly, the observed distribution of A4' is connected to the number density on shells 
of constant z. With a physically plausible assumption regarding the evolution of the source 
distribution function with z, the joint distribution of observed inspiral events as a function 
of A4' and A is sufficient to measure Hq, qq, and otherwise test cosmological models. 



Let MdM-dV be the intrinsic binary coalescence rate (events per time), where M. is 
the intrinsic chirp mass and dV is the local volume element. The SNR p is a function of the 
parameters that describe the binary system and its orientation with respect to the detector. 
We give the form explicitly for a single interferometer in FC (eqs. 3.29, 3.31, and 4.10). The 
quantities dV and H depend on the cosmology and the binary source evolution, and we 
denote the full set of model cosmology parameters as 

The observer counts as detections only those events for which the SNR exceeds a 
threshold po; measuring Al', p, and (for observations in more than one LLI) additional 
parameters C, that are functions of the relative orientations. We wish to compare the 
observed distribution with that implied by a model /x. The expected event rate density 
corresponding to the model is 



3. 



Method 



n 



{M',pX) = / 



dMdV 



dQe M 
An 1 + z 
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5[p - p{dL, M', fiB, fisMM' - (1 + z)M]5[C - CifiB: ns)]. (2) 

The integration over dQ^dQ effectively averages over fiB and hs since, by assumption, J\f 
depends on neither. The total rate n{> po) is h{M', p, C) integrated over all M', C; and 
p > Pq. The probability density for observations exceeding the threshold for a given set of 
cosmological parameters is 

,(M',pX\>Po,f^) = ^^. (3) 

n{> Pq) 

Given a set of N observations we now apply a standard maximum likelihood analysis 
(Eadie et al. 1971). If the confidence volume dM.' dpdVc_ is small compared to the scale on 
which the probability density p varies, then the total likelihood A(/i) oc IljLi (p dM' dp dVc)^ 
The cosmological parameters ^jl that best fit the observed distribution are those that maximize 
A. Regions with A > Aq (for fixed Aq) are confidence volumes for /a. 



4. Application 



As a concrete application, we show how observations in a single LLI can determine 
the parameters of a matter-dominated FRW model. Write the volume element dV as 

f{z] Ho, qo)dzdQ, noting / = f{z; qo){c/ Hof and di = diiz] qo){c/Ho) (where / and di are 
dimensionless and independent of Hq) . 

Assume that Af can be written as Af = N'o{M)h{z){l + z)^ where Af(){M) is the 
intrinsic rate density for coalescences in the local neighborhood. The factor h{z){l + z)^ 
expresses the time-dependence of the intrinsic rate and h{z) encompasses all evolutionary 
effects for a comoving volume [h{0) = 1]. The factorization of J\f implies that the chirp mass 
distribution is independent of the age of the universe, although the coalescence rate may 
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vary. This assumption is justified if the character of the binary formation process (suitably 
averaged over many galaxies) is independent of z. The assumption is very plausible for 
neutron star binaries. Neutron stars form in the collapse of degenerate stellar cores in 
type II supernovae or by accretion induced collapse of white dwarfs. In either case the 
mechanism that triggers the collapse and determines the mass depends largely on the 
equation of state of degenerate matter. The resultant mass distribution should be relatively 
insensitive to the specific evolutionary pathway. The assumption is moderately plausible for 
binaries containing a black hole. Current understanding suggests that the hole's upper mass 
is not tightly constrained by microscopic physical processes and, in fact, observations show 
a range of candidate masses (Casares, Charles and Naylor 1992; Remillard, McClintock 
and Bailyn 1992). The assumption that source evolution is a function of z alone is a 
phenomenological, not a physical, treatment. In contrast, a physical approach would 
introduce additional, non-cosmological timescales. Then h would depend, not only on 
but also go and the timescales for gas cooling, star formation, and so forth. Here, hi^z^ 
serves as a parameterization of the effects of evolution. 

Now assume that neutron star binaries are the predominant sources. Noting that 
the measured range of pulsar masses is small (iO.lM©; Lamb 1991), while the derived 
range of masses in binary X-ray pulsars is somewhat larger (±0.4Mq; Nagase 1989), take 
A/q (TW) = A/q b (TW — M.q)^ i.e., assume all binary systems have the same intrinsic chirp 
mass A^o (we relax this assumption in our more detailed report). Since 2 is a function of 
A^'/A^o, we can express n{z,p) as n(A4', p)/A4o. Using FC (eq. 3.16) for p in eq. (0), we 
find 

dh (z; go) 



(4) 



1 + zf 

where P{x) is the fraction of randomly oriented binaries with B exceeding x [tabulated in 



h{z\ > po) = A^47r (^-|-) f{z; qo)h{z)il + zfP 



FC, where it is called P(G > x)]. The distribution is cutoff as the argument of P increases: 
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P(0) = 1 and P{x > 4) = 0. Note the appearance of the dimensionless factor 

^ = (lOOkms-°Mpc-0 (tI) ■ 

The interferometer's sensitivity is described by F7/3 (for details see FC) which has a nominal 
value of 1 for a single LLI of advanced design and, as used here, -\/2 for the LD composed 
of two such LLI. We refer to the parameters in eq. @ as fiducial parameters. As 5 — 00 
the fraction of observable binaries goes to zero. The total observed rate above threshold is 
eq. @ integrated over all z. 

Assuming no source evolution [h{z) = 1], expand fi{z\ > po) to lowest order in z to 
obtain the Euclidean rate 

n^(>Po)=:^4vr(j-)'(l:^), (6) 

independent of Hq and oc above threshold. We summarize observation rates in terms 
of Q{qo, S, h) = n{> po)/nE{> po). Note lim^^oo Q = I- 



5. Results 



Consider first the characteristic depth of the observed sample. Let Zi/2 (zmax) be the 
redshift at which P = 0.5 (0.0) in eq. (^). To lowest order in redshift, zi/2 = 0.145(9.0/5) 
and Zmax = 0.44(9.0/(5). Note that we have scaled 6 to the fiducial parameters of eq. (|). 
The observed sample is more precisely characterized by the cumulative distribution of 
detected binaries and by its completeness, displayed in figure la for 6 = 4.5, 9.0, 13.5 and 
18.0 (SNR thresholds of 4.0, 8.0, 12.0 and 16.0 for fiducial parameters). Koo and Kron 
(1992, KK) model the distribution of faint galaxies without source evolution but with a low 
Hq. To compare with KK adopt Hq = 50 km/s/Mpc; then the median observed binary has 
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z — 0.10 at po = 8 = 18.0). In KK, the median of the cumulative galaxy distribution 
with 18 < S < 19 occurs at roughly the same redshift. The depth is modest, but (as figure 
la shows) the extent is large. The cumulative distribution for different values of go show 
discernible changes. 

The detection rate n(> po) scales directly with the local coalescence rate. Phinney 
(1991) has discussed the observational and theoretical limits. Based on the several observed 
neutron star binaries, the time until coalescence, and the effective Galactic volumes 
sampled, he gives a "best guess" estimate for the coalescence rate density of 8 x 10^^ Mpc~^ 
yr~^. This is a conservative estimate because it counts only binaries like those observed to 
date; for a number of reasons the true value may exceed the estimate by as much as ~ 800 
[ifo/(100 km/s/Mpc)]^. Scahng to the conservative value, the anticipated detection rate is 

Recall that the function Q is the ratio of the cosmological rate to the Euclidean rate; figure 
lb gives Q for a range of 5 for values of = 0.25, 0.50, and 0.75 assuming no source 
evolution. For the fiducial parameters, Q — 0.715 for go = 0.5, which corresponds to an 
observed rate of h{> 8) = 49. yr""-*^ events. Notice Q is only weakly sensitive to Qq. in 
an analytic expansion to lowest order in z it is independent of go- The effects of source 
evolution may dominate the influence of go on the value of Q. 

Our main interest here is to assess our cosmology test when source evolution takes 
place; consequently we adopt an ad hoc model h{z) = exp[— (a^; + f^z"^)] for z < 10 and 
h{z) — otherwise. This model is useful for both strong and weak evolution. Broadhurst 
et al. (1988) consider a similar form for galaxy luminosity evolution; however, we make 
no attempt to apply their results. In adopting this model, we envision three qualitatively 
different possibilities: (1) The coalescence rate may be simply proportional to the star 
formation rate. Then, based on the faint galaxy studies, there should be little evolution for 



- 10 - 



z < 0.5 (KK). We assume a = /? = for uniform production (no evolution). (2) Neutron 
star binaries may be formed by coUisional encounters in the dense cores of galaxies. The 
luminosity evolution of QSOs is a strong function of z (Boyle et al. 1987, 1988); hence, if 
the coalescence rate is proportional to the QSO luminosity, a substantial increase in the 
production rate with z is expected. We take a = —3/2 and (3 = 1/4 [for which h{z) peaks 
at 2; = 3, where it is /i ~ 9.5] for this case of early production. (3) Binaries may be formed 
by coUisional encounters in globular clusters after core collapse. The time to core collapse is 
of order I/Hq for clusters in the Galaxy. For the case of late production, we take a = and 
(3=1, which cuts off strongly for z > 1. Figure lb illustrates how early and late production 
alters Q in the go = 0.5 model. 

We used a Monte Carlo analysis to calculate the precision with which fj, can 
be determined. We chose "true" model parameters ft and formed a realization of N 
observations. Since we have restricted our analysis to large po the small observational 
uncertainties have not been included in the realization.^ We applied a maximum likelihood 
procedure to this realization to deduce the most likely values jl, and found the difference 
5/x = /i — /i as a function of by repeating the procedure 100 times for each N. For each 
parameter fii, we found the bias {Sfii) and the rms dispersion a^. = (6fii ) . In all cases 
we found the magnitude of the bias to be much less than the dispersion; so, we deal with 
the latter as an indication of the best performance that can be expected for measurement 

of fJ,. 

Our model depends upon five physical quantities: Bubble's constant, the deceleration 
parameter, the chirp mass, a and /5. Denote the "true" values by {Ho,qo,AiQ,a,(3} and 
write the cosmological parameters as {/x} = {ho = Ho/Ho,qo,m = M./J^o,a, P}. (In this 
dimensionless form {/i} = {l,go) Ij^,/?}.) We carried out the Monte Carlo calculations 

^FC show that the fractional errors dAi'/A4' ~ 10~^/p and dp/p ~ 1/p. 
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of the rms dispersion of and go for Einstein-DeSitter cosmologies {i.e., % = 0.5). Our 
results are unchanged by scaling Hq, A4q, po and F7/3 for fixed 6. We assumed uniform a 
priori probabilities for ho G (0.1,2.0) and go ^ (0.1,2.0); the allowed ranges of m, a, and /3 
were large and will be discussed in more detail elsewhere. 

Table |I] summarizes results for and for two values of S, three choices of a and 
P, and several different A^. While the dispersion and bias in ho, go, « and P are all strongly 
correlated, we defer a more detailed discussion to a forthcoming paper. The dispersions 
reported in Table 1 should only be taken as rough approximations: we have not attempted 
to determine the precision to which they are found by our Monte Carlo analysis. It is 
clear, however, that our samples (for fixed 6) extend into the asymptotic regime, where the 
dispersion is oc iV^^/^. 

Note how, at fixed N, the Hubble constant is not better determined by the deeper 
sample (smaller threshold po, which implies smaller 6): thus, sample size is more important 
than sample depth to the accurate determination of ho. In contrast, the uncertainties in 
go are significantly smaller in the deeper sample. More quantitatively, we find that with a 
threshold po ~ 10.7 {6 = 12), LIGO will determine Hq to ~ 10% with ~ 100 observations. 
Similarly, go may be determined to approximately ~ 20% with ~ 3000 observations. 

6. Discussion and conclusions 

Observations of the gravitational radiation from binary inspiral by the LIGO detector 
can be used to measure important cosmological parameters. These measurements rely on 
comparing the observed distribution of inspiral events, as a function of signal strength 
and chirp mass, with model predictions. It is necessary to adopt a model for the intrinsic 
coalescence rate of binaries as a function of the chirp mass and the redshift. We argue 
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that a factorized form is plausible but that other forms are possible. We provide a simple 
demonstration of the power of the test by specializing to the case of a narrow, well-defined 
chirp mass range. We characterize the observational volumes by the sample's median 
redshift and by its completeness. We give the rate of detections as a function of the 
detector threshold signal-to-noise ratio using Phinney's (1991) conservative estimate of the 
local binary coalescence rate. We find that, with a modest number of detections, these 
observations can determine the Hubble constant to 10%. The conservative rate estimate 
imphes that the observations would take ~ 4 years. We address the optimal choice of the 
SNR threshold, measurement errors at small p and a range of binary masses in a paper in 
preparation. 

We thank Drs. Bildsten, Cutler, Flanagan, Loredo, Thorne, Wasserman, and Whitcomb 
for discussions. LSF acknowledges the Sloan Foundation and NASA NAGW-2936 
(Northwestern); DFC the PYI program, NSF AST-8657467 and NASA NAGW-2224 
(Cornell). 
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FIGURES 

Figure la. The observational volume accessible to the LD is characterized by the 
cumulative distribution of detected hinaxies S{z,pq) = Jq dz'h{z'\ > po)/h{> po) (the rising 
set of curves) and by the sample's completeness T{z, pq) = Jq dz'h{z'\ > po)/ lo dz'h{z'\ > 0) 
(the falling curves) for the case of no source evolution. The widely spaced families of 
curves represent 6 = 4.5, 9.0, 13.5 and 18.0 (eq. ^ from right to left; these correspond to 
thresholds of po = 4, 8, 12 and 16 for fiducial parameters. Each family consists of three 
closely spaced curves corresponding to different values of go (dashed 0.25; solid 0.5; dotted 
0.75). The curves describing the sample's completeness terminate at the maximum z at 
which an inspiraling binary can be observed. 

Figure lb. The ratio of the cosmological to the Euclidean rate of detections, Q, is 
given as a function of S (eq. For the case of no source evolution, three different values 
of go are displayed, as in Figure la; these correspond to the three closely spaced lines. 
Two additional models discussed in the text are shown for the go = 0.5 case. The long 
dashes illustrate a coalescence rate lower in the past than in the present. The chain dashes 
illustrate one with a higher rate in the past. 
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Table 1: The results of applying the maximum-likelihood analysis to finding ho, q^, ttiq, a 
and (3. 

5 N No Evolution Early Production Late Production 
(a = /3 = 0) (a = 0, /3 = 1) (a = -1.5, /3 = 0.25) 

f^ho CTqo ^ho CTqo ^/lo <^?o 



30 


0. 


,144 


0.844 


0. 


,143 


0.838 


0, 


.130 


0.881 


100 


0. 


,099 


0.681 


0. 


,093 


0.674 


0, 


.104 


0.666 


300 


0, 


,071 


0.479 


0, 


,063 


0.425 


0, 


,075 


0.425 


1000 


0. 


,035 


0.197 


0. 


,039 


0.226 


0, 


.036 


0.193 


3000 


0. 


,019 


0.110 


0. 


,017 


0.102 


0, 


.020 


0.110 


10000 








0. 


,009 


0.055 


0, 


.010 


0.054 


30 


0. 


,097 


0.890 


0. 


,112 


0.865 


0, 


.109 


0.896 


100 


0. 


,075 


0.780 


0. 


,076 


0.742 


0, 


.067 


0.696 


300 


0. 


,054 


0.590 


0. 


,049 


0.515 


0, 


.054 


0.558 


1000 


0. 


,029 


0.311 


0, 


,026 


0.255 


0, 


.027 


0.273 


3000 


0. 


,013 


0.141 


0. 


,015 


0.157 









